Fluctuation relation for a Levy particle 



H. Touchette^B and E. G. D. Cohen 2 

1 School of Mathematical Sciences, Queen Mary, University of London, London El ^AW, UK 
2 The Rockefeller University, 1230 York Avenue, New York, New York 10021, USA 

(Dated: February 6, 2008) 

We study the work fluctuations of a particle subjected to a deterministic drag force plus a random 
forcing whose statistics is of the Levy type. In the stationary regime, the probability density of the 
work is found to have "fat" power-law tails which assign a relatively high probability to large fluc- 
tuations compared with the case where the random forcing is Gaussian. These tails lead to a strong 
violation of existing fluctuation theorems, as the ratio of the probabilities of positive and negative 
work fluctuations of equal magnitude behaves in a non-monotonic way. Possible experiments that 
could probe these features are proposed. 
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The study of fluctuations in systems driven in nonequi- 
librium steady states has concentrated lately on a special 
symmetry property referred to as the fluctuation theorem 
or fluctuation relation. The first studies of this symme- 
try focused on the entropy production of nonequilibrium 
systems 0, but it was soon realized that the same 
symmetry holds for other quantities of interest, such as 
the work W T performed on a nonequilibrium system dur- 
ing a time r [3| . In the context of this quantity, we define 
a fluctuation relation Q as follows. Assuming that W T is 
extensive in r, we say that the probability density P(W T ) 
of W T satisfies a fluctuation relation or, to be more pre- 
cise, a conventional fluctuation relation if 

P(W T = wt) _ 

P(W t = -wt) ' {) 

where c is a constant which depends neither on t nor w. 
We also say that P(W T ) satisfies a conventional fluctua- 
tion relation if the above equality holds, when properly 
scaled (see below), in the limit r — > oo. In that case, we 
speak of a stationary conventional fluctuation relation. 

One important problem about fluctuation relations, 
which is as yet unresolved, is to determine the class of 
nonequilibrium systems or, more precisely, the class of 
nonequilibrium observables [f| whose fluctuations satisfy 
a conventional fluctuation relation. Our understanding 
of this issue at this point is that this relation holds for 
two general classes of systems: (i) finite, deterministic 
systems having a sufficiently chaotic dynamics @, [(| , and 
(ii) finite, stochastic systems whose evolution is a Markov 
process [7|- There are subtle points to be taken into ac- 
count, however. In some cases, boundary conditions or 
special forms of large deviations may restrict the range 
of validity of conventional fluctuation relations, and in 
these cases, corrections or extensions of the conventional 
fluctuation relation have been proposed [1, Q . 

Our goal in this paper is to expand this picture of fluc- 
tuation relations by studying a model of a nonequilibrium 
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system whose work fluctuations neither obey the conven- 
tional fluctuation relation nor the extended fluctuation 
relation of van Zon and Cohen Q . Based on this model, 
we propose a novel type of fluctuation relation, and com- 
pare it, from the general point of view of large deviation 
theory, with the conventional fluctuation relation defined 
in ([!]). In the end, we also discuss two experiments for 
which "anomalous" fluctuations similar to the ones found 
for our model could be observed. The first experiment is 
related to dry friction, while the second draws on recent 
experiments on granular gases. 

The specific model that we study is a variation of the 
Langevin model studied recently by van Zon and Cohen 
Q, which describes the motion of a particle subjected to 
three forces: a deterministic force derived from a particle- 
confining harmonic potential, a friction force, and a ran- 
dom force, i.e., a noise. The combined effect of these 
forces on the particle is described, in the overdamped 
limit, by the Langevin equation 

i(t) = -7b(*) "**(*)] +£(*), (2) 

where x(t) denotes the position of the particle at time t, 
x* (t) the position of the center of the harmonic potential 
at t, 7 a constant related to the ratio of the strengths 
of the potential and the friction, and £(t) the random 
forcing. For simplicity, we assume that the center of the 
potential moves at constant speed v*, i.e., x*(t) = v*t. 
Finally, and this is where we depart from the model of 0] , 
we assume that £(t) represents a Levy white noise, i.e., 
an uncorrelated noise obeying a Levy statistics [l(| ■ Such 
a noise process is entirely defined by its characteristic 
function 

G ( [k] = (e* k <) = J P[f] e'XT ( 3 ) 

which, in the case of s ym metric Levy white noise, is taken 
to have the form [TolTrH 

G 6 [k] = exp ^-6^ \k(t)\ a dtj (4) 

with b > and < a < 2. Gaussian white noise is a spe- 
cial case of the noise process, obtained by choosing a = 2. 
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For this case, which is the typical case considered by van 
Zon and Cohen, (£(*)) = and (£(*)£(*')} = 2 ^(* ~ *')■ 

The introduction of Levy noise with a € (0, 2) in 
Langevin equations is often perceived as being physically 
unreasonable because it causes x(t) to have an infinite 
variance, which implies that there can be no fluctuation- 
dissipation relation connecting the strength of the friction 
to the variance of the noise [10]. This criticism indeed 
applies if we regard the noise as being internal, i.e., as 
arising, like the friction, from the environment surround- 
ing the particle. However, if we regard the noise as be- 
ing external, i.e., as being an exogenous driving force, 
then there is no problem in assuming that is Levy- 
distributed. In that case, the friction and the external 
noise are decoupled, and one is free to adopt any model 
of noise. This point will be clarified again when we dis- 
cuss the experiments that we have alluded to above. 

For now we leave these considerations aside, and at- 
tempt to calculate the probability density P(W T ) of the 
work W T performed on the Levy particle by the harmonic 
drag force during a time r 



W T = -v* / [x(t) -x*(t)]dt. 
Jo 



(5) 



To carry out this calculation, it is convenient to study 
the dynamics of x{t) in a comoving frame with the 
change of coordinate y(t) — x(t) — x*(t), so that W T — 
—v* J Q y(t)dt. The dynamical equation for y(t) reads 



m = -iy(t) + r,(t), 



(6) 



where r)(t) = £(t) — v* represents a white noise with mean 
or "drift" — v* . 

At this point, we proceed to calculate P(W T ) by cal- 
culating its characteristic function, defined as 

poo 

G W M = (e iqWr ) = / P{W T = w) e iqw dw. (7) 



A key mathematical observation here is that this charac- 
teristic function can be obtained from the characteristic 
function G y [k] of the process y(t) by using k(t) = —v*q 
for t € [0,t] and k(t) = otherwise. With this choice of 
k(t), we indeed obtain from Eqs. (TJJ and §5§ 

G v [k] = /cxp (i [ k(t)y(t)dt)\ = G Wr (q). (8) 



In order to find G y [k], we next solve the Caceres-Budini 
formula associated with y(t) [l2T ]: 

/oo 
e^ l -*h{t)dt. (9) 

In applying this formula to our model, we assume that 
j/(0) = x(0) = 0. We also adopt at this point dimension- 
less units by setting 7 = 1 and 6 = 1. Noticing that 



G v [r] = cxp ( -iv* j r(l)dl ) G 6 [r] , (10) 
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FIG. 1: (Color online) (a) Plot of P(W T = w) for various 
values of a (M = V = 1). (b) Log- Log plot of P(W r = w). 



we can write 

/ poo poo \ 

G WT (q)=expl-iv* J r(l)dl-j \r(l)\ a dlj . (11) 

By performing the integrals involved in this expression, 
we arrive at 



G W M) = e lMq - v ^ c 



where 



M = u(t + er T - 1) 



V 



/Q /2 



I — T 



ir di 



(12) 

(13) 
(14) 



having defined v = (v*) . 

This characteristic function is the main mathematical 
result of this paper, from which P(W T ) is obtained by an 
inverse Fourier transform. In the present case, there is no 
need to explicitly compute the inverse Fourier transform 
because the form of Gw T {o) shown in (fTTj) already tells 
us that P(W T ) is a symmetric Levy distribution having 
the same index a as the noise £(t) [ll|. The parameter 
M is the value at which P(W T ) is centered, while V is a 
measure of the width of P(W T ); see Fig. [1] The index a, 
finally, determines how the tails of P(W T = w) decay as 
\w\ — > 00. It is at this point that we encounter two very 
different types of fluctuations. 

A first type arises from Gaussian white noise, which 
is obtained again with a = 2. For this case, and for the 
particular initial condition that we consider here, namely 
x(0) — 0, P(W T ) has an exact Gaussian form 



P(W T =w) = 



1 



V27TCT 2 



-(w-M) 2 /(2a 2 ) 



(15) 



with mean M given by (|13[) and variance a 1 given by 

a 2 = 2V = v[2r- 3 + 4e- T - e~ 2r ]. (16) 

This form of P(W T ) obeys a stationary conventional fluc- 
tuation relation. Indeed, since M ~ vr and cr 2 ~ Ivt ~ 
2M in the limit r — > 00, we have 



-(P) 



P(W T = pvr) 
P(W T = -pur) 



(17) 
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FIG. 2: (a) Comparison of g T (jp) and (b) \ng T (p) for the 
Cauchy (solid line) and Gaussian (dashed line) cases (a = 1, 
1/ = 1, and t — 1 for the Gaussian case). 



with sub-exponential corrections in r. This result con- 
firms what van Zon and Cohen have reported in Q using 
a different calculation technique. Here, p is the scaled 
(intensive) value of W T defined as p — W T /(W T ), with 
the brackets denoting an ensemble avarage. 

An entirely different case of fluctuations arises when 
a < 2, i.e., when a £ (0,2). In this case, P(W T ) is not 
Gaussian an ym ore. Rather, it is a so-called strict Levy 
distribution [ll| having, as is well known, power-law tails 
of the form 



P(W T = W) 



(18) 



as an ±oo; see Fig. QJb). Because of these power-law 
tails, it is clear that P(W T ) does not satisfy a fluctuation 
relation of the form defined by Eq. fl}, even in the long 
time (stationary) regime. To be sure, consider the case 
a = 1, for which we have the following exact solution: 



1 



P{W T = w) = - 



V 



n(w- M) 2 + V 2 ' 



(19) 



with M given, as before, by Eq. (|13[) and V — zW 2 (r — 
I + e~ T ) ^ v x I 2 t. This form of distribution is known as 
a Cauchy or Lorentz distribution, and behaves like |w|~ 2 
for large work values. Using this explicit form of P(W T ), 
we obtain in the limit r — > oo 



9r(p) 



v(p+l) 2 + l 

Kp-i) 2 + i' 



(20) 



As shown in Fig. [21 this form of g T (p) is not a mono- 
tonic function of p, contrary to its Gaussian counter- 
part. In the Cauchy case, we further have that g T {p) does 
not asymptotically depend on time, which implies that 
measuring W T over longer and longer integration times 
r leads to no appreciable difference in the behavior of 
the fluctuations of the scaled variable p. In the Gaussian 
case, g T {p) does depend on time, with the consequence 
that negative fluctuations of the work get exponentially 
suppressed as r — > oo. 

The asymptotic time independence of g T {p) for Cauchy 
white noise is special, and does not carry over to other 
values of a in the range (0,2). However, since all strict 



Levy distributions have power-law tails, P(W T ) cannot 
satisfy a conventional fluctuation relation for any a £ 
(0,2). Indeed, because of (fl8l) . we have 



lim g T {p) = 1 

p — '±00 



(21) 



for all a £ (0, 2) (Fig. [2|). In words this means that large 
negative fluctuations of W T are just as likely to happen as 
large positive fluctuations of equal magnitude. For a = 2, 
by comparison, we have g T (p) — > oo when p — ► oo and 
g T (j>) — > when p — > — oo, which implies again that pos- 
itive fluctuations of W T are exponentially more probable 
than negative fluctuations. This difference in behavior 
clearly shows that the fluctuations of the work performed 
on a particle driven by strict Levy white noise satisfy an 
entirely different kind of fluctuation relation when com- 
pared with its Gaussian noise counterpart. Following a 
terminology used in studies of Levy noise [TlT |. we call 
this new fluctuation relation, satisfying both (fT8")) and 
pip , an anomalous fluctuation relation. 

An important property of the density P(W T ) brought 
about by its power-law tails is that it does not have the 
form of a large deviation probability. By this we mean 
that P(W T ) does not have the asymptotic form 



P(W T = wt) ps e 



-tI(w) 



(22) 



in the limit r — > oo, where I(w) is some function which 
does not depend on r fl3l |. This observation is crucial 
because it explains why P(W T ) fails to satisfy a normal 
(i.e., conventional or extended) fluctuation relation even 
if the Langevin model studied here is Markovian. The fact 
is that the Markov property is not a sufficient condition 
for having a normal fluctuation relation, as is demon- 
strated by our model. What is essential is that the ob- 
servable of interest be governed by a probability density 
having a large deviation form. Note indeed that in order 
for P(W T )/P(—W T ) to be exponential in r, P{W T ) must 
itself be exponential in r, as in (f2"2"]) . Here, P(W T ) has 
a large deviation form only for a = 2, which is why this 
value gives rise to a normal fluctuation relation, while all 
the other values of a give rise to an anomalous fluctua- 
tion relation. A more detailed explanation of this point 
together with a complete analysis of P(W T ) for values of 
a other than 1 and 2 will be given in another paper [l4[ . 

To complement the theory presented so far, we now 
suggest an experiment that could be used to observe the 
anomalous fluctuation relation displayed by our model. 
As mentioned before, the type of physical system that 
we must look for is one in which the noise is introduced 
externally and can be tailored to mimic Levy white noise. 
One such system is provided by an experiment recently 
reported by Buguin et al. [HI], in which a solid object 
placed on a solid plane is set in motion by rapidly vi- 
brating the plane horizontally. In [l6| , de Gennes consid- 
ered the possibility of vibrating the plane using Gaussian 
white noise, and provided various theoretical predictions 
for this case, based on a Langevin equation containing 
a dry friction term and a noise term, uncorrelated with 
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the friction term, representing the random motion of the 
plane. Repeating the experiment of Buguin et al. with 
Gaussian white noise, one could check whether a nor- 
mal fluctuation relation for the work W T performed by 
the plane on the moving object during a time r holds. 
Then, using Levy white noise, one could check whether 
an anomalous fluctuation relation holds for W T [l7j . 

A different experimental investigation of anomalous 
fluctuations could be based on a gas of inelastically- 
colliding particles, maintained in steady states by rapidly 
vibrating its container. Several studies of this so-called 
"granular gas" have appeared recently (see, e.g., [H|,[l||), 
and focus on the fluctuations of the power injected by the 
vibrations. In this context, two types of vibrations have 
been considered, namely periodic forcing and Gaussian 
white noise, and for these it appears that normal fluc- 
tuation relations hold in the short time regime where r 
is comparable with the particles' mean free time [l8|. 
In view of our model, one could change the vibrations 
from Gaussian to Levy white noise, and see whether an 
anomalous fluctuation relation similar to the one found 
here holds for the short time fluctuations of the injected 
power. One could also check that the distribution of the 
injected power does not have the form of a large devia- 
tion probability when strict Levy noise is used. This can 
be done in simulations, and, at least in principle, in real 
experiments. 

Note that Levy noise, like Gaussian noise, is a proper 
form of noise, in the sense that its realizations are always 
finite [11] . Therefore, forcing an object with a Levy noise 



having an infinite variance, or even an infinite mean, does 
not imply that we supply that object with an infinite 
amount of energy. In particular, Levy noise cannot sup- 
ply an infinite amount of heat to a granular gas in a finite 
amount of time. 

To summarize, we have shown that the work fluctua- 
tions of a nonequilibrium system can show different be- 
haviors depending on the form of noise perturbing the 
system. This was illustrated with a Langevin equation 
perturbed by Levy white noise. For this model, we have 
shown that the probability density of the work fluctua- 
tions does not have the form of a large deviation probabil- 
ity, since it has power-law tails. These fluctuations lead 
therefore to a novel type of fluctuation relation, which 
differs from the conventional and extended fluctuation 
relations discussed up until now. Other noises behaving 
asymptotically like Levy white noise should give rise to 
similar results, since Levy distributions are stable [Tl| . In 
this sense, the results obtained here should be represen- 
tative of a large class of fluctuation relations associated 
with "power-law" noises having infinite variances. 
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